Surface soft phonon and the \/3 x a/3 ^3x3 phase transition in Sn/Ge(lll) and 

Sn/Si(lll). 



O 

o 

(N 



c3 



-a 
c 

o 
o 



> 

(N 



O 

o 

-I— > 
c3 



i 

c 

o 
o 



X 



Ruben Perez, Jose Ortega and Fernando Flores 
Departamento de Fisica Teorica de la Materia Condensada, Universidad Autonoma de Madrid, E-28049 Madrid, 
(manuscript LY7332, final version accepted in Physical Review Letters) 

Density Functional Theory (DFT) calculations show that the reversible Sn/Ge(lll) a/3 x a/3 ^3x3 
phase transition can be described in terms of a surface soft phonon. The isovalent Sn/Si(lll) case 
does not display this transition since the a/3 x a/3 phase is the stable structure at low temperature, 
although it presents a partial softening of the 3x3 surface phonon. The rather flat energy surfaces 
for the atomic motion associated with this phonon mode in both cases explain the experimental 
similarities found at room temperature between these systems. The driving force underlying the 
a/3 x a/3 <-» 3 x 3 phase transition is shown to be associated with the electronic energy gain due to 
the Sn dangling bond rehybridization. 



PACS numbers: 73.20.-r, 73. 20. At, 71.15.Nc, 71.15.Ap 



A soft phonon is a vibrational mode of a crystalline 
material whose frequency decreases as T falls, eventually 
reaching zero. At this point the crystal is unstable in 
relation to the corresponding atomic displacements and 
undergoes a transition to a lower symmetry phase. Typ- 
ical examples are the high-low transition in quartz, the 
ferroelectric transition in BaTiO^ and the ferroelastic 
transition in SrTiO?, Q|. On semiconductor surfaces, 
different mechanisms, like Peierls transitions or Charge 
Density Waves, have been proposed as the origin of two 
dimensional phase transitions. However, no soft phonon 
transition has been unambiguously identified in these sys- 
tems so far. 

In this letter we consider the Sn/Ge(lll) a/3 x a/3 
3x3 phase transition, which is a prototypical example of 
a T-induced reversible change of symmetry at a semicon- 
ductor surface. The driving force underlying this phase 
transition has been under intense debate since its dis- 
covery on the Pb/Ge(lll) system [0. Different physical 
mechanisms have been proposed to play an important 
role in this transition: surface Fermi wavevector nest- 
ing or electron correlations, leading to the formation of 
a Charge Density Wave |^,|| at low T; dynamical fluc- 
tuations of an underlying 3x3 structure, giving rise to 
the observation, on average, of a a/3 x a/3 symmetry at 
room T B ; and the interaction between the 3x3 period- 
icity and Ge-substitutional defects that act as nucleation 
centers 0]. Recently, several groups have also con- 
sidered the isovalent Sn/Si(lll) system and have found 
strong similarities and some striking differences with the 
Ge case. The electronic structure is very similar in both 
cases, with two major components in the Sn Ad core level 
spectra and two surface bands close to the Fermi level 
M, associated with two different Sn-dangling-bonds. At 
variance with the Sn/Ge(lll), where a transition to a 
3x3 pattern is clearly observed at low T, the Sn/Si(lll) 
surface shows a a/3 x a/3 pattern both in LEED and STM, 
for temperatures as low as 70 K S, even in the presence 
of a significant density of Si-substitutional defects 0. 



Moreover, the two surface bands mentioned above are not 
as clearly resolved at low T as they are in the Sn/Ge(lll) 
case ||. By analyzing theoretically the Sn/Ge and Sn/Si 
systems, we conclude that the a/3 x a/3 3 x 3 tran- 
sition is due to a surface soft phonon. The frequency 
of this mode in Sn/Ge(lll) goes to zero for a fc-vector 
corresponding to a reciprocal lattice vector of the 3x3 
periodicity. For Sn/Si(lll) we find that the a/3 x a/3 
is the stable structure at low T, and therefore this sys- 
tem should not display the a/3 x a/3 -^3x3 transition. 
The comparison of the Sn/Ge and Sn/Si systems with 
their corresponding H-covered cases shows that the driv- 
ing force for this transition is the electronic energy gain 
associated with the surface band splitting induced by the 
3x3 distortion. 

The atomic displacements associated with the a/3 x 
a/3 3 x 3 transition in Sn/Ge(lll) provide the clue 
for the soft mode responsible for this transition. The sig- 
nature of the 3x3 structure is the upward displacement 
of one of the three Sn-adatoms, and the corresponding 
downward movement of the other two JlOf . Figure 1 and 
Table I show how these upward and downward move- 
ments are interconnected: as one of the Sn atoms (Sni) 
moves upwards, the three Ge nearest neighbours (A\) fol- 
low this motion by moving towards it (with both upward 
and in-plane displacements) . These Ge displacements 
force corresponding in-plane movements of the Ge (la- 
belled D2) towards the position of Sni. Since the Ge-D2 
atoms are also bonded to first-layer Ge-atoms (Bi,Ci) 
linked with the two other Sn atoms (SV12 and Sn$), these 
are forced to move downwards. 

Does the Sn/Si(lll) surface present a similar 3x3 
structure at low T?. We have studied this possibility 
by means of DFT calculations JO]. In this analysis, we 
start from a a/3 x a/3 structure (in a 3 x 3 unit cell), 
where the 3 Sn atoms are equivalent, we select one of 
them (Sni) and force it to move in the direction perpen- 
dicular to the surface. For each of these displacements, 
the other Sn atoms (Sn% and SV13) and all the semicon- 
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ductor atoms are allowed to relax (up to the fifth layer) 
to their zero force positions under the constraint of the 
Sni displacement. These calculations were performed for 
both the Sn/Si and Sn/Ge systems (see Figure 2). The 
energy vs. Sni displacement is very flat for Sn/Ge(lll); 
the ground state corresponds to a 3 x 3 geometry with a 
Sni displacement of ~ 0.18 A, as found in previous cal- 
culations pX|Jl2| and in agreement with the experimental 
evidence pJLp'j]. The Sn/Si(lll) surface, however, does 
not present any stable 3x3 distortion, i.e., the v3 x \/3 
is the stable structure at low T. 

The inset of Fig. 2 shows another interesting finding: 
Sn,2 and Sn$ move together in the opposite direction to 
the constrained motion of the Sni atom, for both the 
Sn/Ge(lll) and Sn/Si(lll) surfaces. These results sug- 
gest a correlated up/down motion of the Sn atoms ||] 
that, at a sufficiently low temperature, is frozen into the 
3x3 structure in the case of Ge (but not in the case of 
Si). As reflected in Table I, the displacements associated 
with this motion are essentially localized in the tetrahc- 
dra formed by a Sn atom and the three Ge of the first 
layer bonded to it. In order to analyze the correspond- 
ing normal modes associated with the up/down motion 
of the tetrahedra, we have considered a force constant 
model which includes effective interactions between first 
nearest neighbour tetrahedra: 

Fj = -azj + (zj - Zj) (1) 

j,nn 

where Fi, the force acting on each tetrahedron, is propor- 
tional to its center-of-mass displacement (with respect 
to the V3 x V3 structure), z%, and the differences in dis- 
placement between nearest neighbours, Zj — z%. 

The force constants a and f3 can be determined from 
the calculations presented above. These calculations pro- 
vide the accurate information required for this purpose, 
since they are focused precisely on the deformations we 
are analyzing here. Imposing that Fi and F3 are zero 
in the above equations yields a relation between the 
force constants and the displacements z\ and z%{= Z3): 
3/3/a = — Z2I (zi ~ zi). For small displacements, (—^2) is 
around 0.5zi for Ge, while it is only 0.35 zi in Si. There- 
fore, we obtain 3/3/a = —1/3 for Ge and 3/3/a = —0.23 
for Si. 

The phonon dispersion relation of this mode is easily 
found, using eqn. |l|, to be: 

Mu 2 = {a + QP)-f3^ cos(k ■ Rj) (2) 
i 

where k is the momentum parallel to the surface, Rj the 
coordinates of the six first nearest neighbours, and M the 
tetrahedron mass. In this approach, the atoms outside 
the tetrahedra are treated adiabatically, i.e. we neglect 
their mass and assume them to follow "instantaneously" 



the displacement imposed by the motion of the tetrahe- 
dra. This means that the forces between, say, the Ge-D2 
atoms (see fig. Q) and the atoms in the tetrahedra are 
automatically included in the force constants of eqn. [l| 
and in the phonon dispersion relation of eqn. || We have 
estimated the error introduced by this adiabatic approx- 
imation in the phonon frequencies by comparing the ki- 
netic energies of the atoms whose inertia is neglected in 
our approach with that of the atoms in the tetrahedra. 
This yields errors of 7% and 16 % for the Si and Ge cases 
respectively. These values set the order of magnitude 
of the accuracy that our model presents for the phonon 
frequencies. 

Figures [|a and [|b show the phonon dispersion curves 
for Ge and Si, compared with the projected phonon bulk 
band structure. Notice that we find a zero-frequency 
mode at the K point for Ge. The corresponding k -vector 
defines the new 3x3 periodicity associated with this soft 
mode. For Si, the phonon dispersion relation shows only 
a minimum at the K point, with a phonon energy of 5.5 
meV. These results are in agreement with our previous 
comments on the stability of the Sn/Si(lll)-v / 3 x \/3 
structure, as shown in Figure ^. 

The phonon branches shown in Fig. || provide an expla- 
nation for the main similarities and differences between 
the Si and Ge cases. At low T, while for Ge we can 
expect a phase transition to the 3x3 structure, there is 
no stable 3x3 phase for Si. On the other hand, at high 
T (say, for ksT >> 6 meV, room T is a high T limit) 
we can expect Ge and Si to be alike. This can be un- 
derstood by considering the energy curves in Fig. ^: at 
room T the Sn atoms vibrate with a large amplitude in 
both surfaces (with displacements between -0.2 and 0.3 
A for Ge and -0.2 and 0.2 A for Si). In this vibration, Sn 
atoms are correlated with their Sn-nearest-neighbours. 
This is true even in the Si case, because although the 
3x3 structure is not an stable minimum, due to the fact 
that the vibration associated with the 3x3 structure is 
a minimum of the phonon dispersion, the atomic motion 
can be expected to be dominated by this mode. For the 
Sn/Ge(lll)— v3 x y/3 phase, anharmonic effects renor- 
malize the frequency of the soft-phonon mode to a non- 
zero value [ [l5|| , as shown schematically in Figure 3. 

The atoms involved in this vibration spend most of the 
time around the turning points of their classical trajec- 
tories, where the atom velocity goes to zero. At these 
turning points, the electron occupation of the Sn dan- 
gling bonds is different: atoms in "up" positions have a 
fully occupied dangling bond, while those at "down" po- 
sitions present only a partial occupation of their dangling 
bonds jl(| . This explains the double peak observed in the 
XPS Sn Ad core level spectra of both systems, as origi- 
nally proposed in the dynamical fluctuations model Q to 
explain the Sn/Ge(lll) 3x3^ \/3 x \/3 transition. In 
this model, Sn atoms display correlated up/down vibra- 
tions which present, at RT, large amplitudes, but keep 
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memory of the underlying 3x3 phase. The soft phonon 
shown in fig. |^ is responsible for this atomic motion and 
provides the physical mechanism for the dynamical fluc- 
tuation model. 

We have also addresed the problem of which driving 
force is softening the surface mode associated with the 
3x3 structure. As Si and Ge are isovalent atoms and 
the electronic structure of both interfaces is very similar, 
one might expect the difference in the relative stability 
of the two structures to be related to the different elastic 
properties of the two materials. In order to single out 
the elastic contribution from the one associated with the 
Sn-dangling-bond rehybridization fl6| and correspond- 
ing surface-band splitting fib]] , we have studied the sur- 
face phonons of these surfaces with hydrogen atoms satu- 
rating the Sn-dangling-bonds. As shown in Figs. ||a and 
||b, the dispersion curves for the H-saturated surfaces 
are rather different from the Sn-surfaces; in particular 
they display no "softening" of the mode at K . This 
different behaviour is related to the Sn-dangling-bonds 
rehybridization. 

In order to quantify these effects, we have considered 
the relaxed Sn-surfaces discussed above (see Fig. 2) 
for two different Sni displacements, namely, 0.0 A(the 
\/3 x V3 geometry) and 0.18 A(the "3x3" geometry). 
The elastic contribution is calculated, in each case, as 
the difference in energy between the two geometries 
with H-atoms saturating the Sn-dangling-bonds. Sub- 
stracting this clastic contribution from the total energy 
difference for the Sn-surfaces, the energy contribution 
due to the Sn-dangling-bonds rehybridization is ob- 
tained. These calculations provide the following results: 
EeiasticiGe) = 156 meV and E e i ectronlc {Ge) = -158 
meV; E e i asUc (Si) = 131 meV and E e i ectronic (Si) = -112 
meV [0. Therefore, the main difference between the 
Sn/Ge and Sn/Si cases is the different hybridization be- 
tween Sn and the surface atoms. A comparison of the 
surface bands for both cases confirms this result: the Ge 
case presents in the 3x3 geometry a larger band splitting, 
152 meV as compared with 93 meV for Si. This effect can 
be related to the different sizes of Si and Ge (r| . When 
Sn is closer in size to the semiconductor atom the surface 
rehybridization is stronger and the surface gains more en- 
ergy by means of the surface band splitting caused by the 
3x3 distortion. 

In conclusion, we have shown via state-of-the-art DFT 
calculations that the reversible V3 x \/3 ^3x3 phase 
transition in Sn/Ge(lll) is associated with a surface soft 
mode. The 3x3 structure is not stable in the Sn/Si(lll) 
case, although it is a minimum in the phonon disper- 
sion, and consequently this surface should not display 
the -\/3 x -\/3 ^> 3 x 3 transition. Comparing the Sn/Ge 
and Sn/Si cases, and considering the corresponding H- 
covered surfaces, we have found that the driving force 
underlying the \/3 x V3 «3x3 phase transition lies in 
the electronic energy gain due to the Sn dangling bond 



rehybridization. 

This work has been partly funded by the Spanish CI- 
CYT under contract No.PB-97-0028. We acknowledge 
helpful discussions with V.R. Velasco and J.M. Perez- 
Mato. 



[1] Andrew Putnis, Introduction to Mineral Sciences (Cam- 
bridge University Press, Cambridge, 1992). 
[2] J.M. Carpinelli et al., Nature 381, 398 (1996). 
[3] J.M. Carpinelli et al., Phys. Rev. Lett. 79, 2859 (1997). 
[4] J. Avila et al., Phys. Rev. Lett. 82, 442 (1999). 
[5] H.H. Weitering et al., Science 285, 2107 (2000). 
[6] R.I.G. Uhrberg et al., Phys. Rev. B 62, 1 (2000). 
[7] L. Ottaviano, M. Crivellari, L. Lozzi, S. Santucci, 

Surf.Sci. 445, L41 (2000) 
[8] T. Yamanaka and S. Ino, Phys. Rev. B 61, R5074 (2000). 
[9] A. Charrier et al. to be published. 
[10] J. Ortega, R. Perez and F. Flores, J. Phys.: Condens. 

Matter 12, L21 (2000). 
[11] L.J. Clarke, I. Stich and M.C. Payne, Comp. Phys. 

Comm. 72, 14 (1992). 
[12] S. de Gironcoli et al., Surf.Sci. 454-456, 172 (2000) 
[13] O. Bunk et al., Phys. Rev. Lett. 83, 2226 (1999). 
[14] Jiandi Zhang et al., Phys. Rev. B 60, 2860 (1999). 
[15] D. Vanderbilt, Current Opinions in Solid State and Ma- 
terial Science 2, 701 (1997). 
[16] G. Ballabio et al., Phys. Rev. B 61, 13345 (2000). 
[17] E elasU c(Si) is smaller than E e i a3tic (Ge) due to the 
smaller displacements of the Si-atoms, for the same Sni 
displacement. 



TABLE I. Atomic displacements (in A) associated 
with the Sn/Ge(lll)-(3x3) surface measured w.r.t. the 
V3 x \/3-geometry. 





Sni 


Sn 2 


Sn 3 


Ai 




Bi 


Ci 


Az 
A _L 


0.18 
0.00 


-0.08 
0.00 


-0.08 
0.00 


0.08 
-0.04 




0.04 
0.02 


-0.03 
0.02 








A 2 


B 2 




c 2 




D 2 


Az 
A _L 




0.02 
0.00 


0.00 
0.00 




0.00 
0.00 




0.00 
0.04 
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FIG. 1. Ball-and-stick model of the Sn/Ge(lll)-(3x3) 
reconstruction. The arrows show the direction of the atomic 
displacements (up/downwards on top, in-plane on bottom fig- 
ure) w.r.t. the y/S x \/3-geometry. Notice that the Ge atoms 
move in-plane towards the position of the upper (Sni) and 
away of the two other Sn atoms. The atoms D2 connect the 
Ge atoms which are nearest neighbours of the different Sn 
atoms. 
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FIG. 2. Total energy as a function of the Sni for the 
Sn/Ge(lll) and Sn/Si(lll) surfaces. The energy of the 

structure is taken as the reference. Inset: Displace- 
ments (in A) of the Sn2 , S113 atoms as a function of the 
Sni displacement. All the displacements are referred to the 

structure. 
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FIG. 3. Phonon dispersion curves for the Sn/Ge(lll) 
(full line, top) and Sn/Si(lll)- x ^3 (full line, bot- 
tom) surfaces. Dashed lines correspond to the H-saturated 
Sn/Ge(Si)(lll)-v / 3x \/3 surfaces. Shaded areas represent the 
projection of the corresponding phonon bulk band structure. 
The phonon branches are plotted along the FM direction in 
the first Brillouin Zone (BZ) of the ideal (111) surface (see in- 
set). The inner hexagon corresponds to the BZ of the y/3 x \/3 
surface. K' defines the new 3x3 periodicity associated with 
the soft phonon in Ge. At room T, the frequency of the soft 
mode should be renormalized to a non-zero value, as shown 
schematically by the dotted line. 
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